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Abstract 



The adhesion between a polymer melt and substrate is studied in the 
presence of chemically attached chains on the substrate surface. Extensive 
molecular dynamics simulations have been carried out to study the effect of 
temperature, tethered chain areal density (S), tethered chain length (A't), 
chain bending energy (kg) and tensile pull velocity (v) on the adhesive failure 
mechanisms of pullout and/or scission of the tethered chains. We observe a 
crossover from pure chain pullout to chain scission as A'^t is increased. Below 
the glass transition, the value of Nt for which this crossover begins approaches 
the bulk entanglement length N^- For the values of Nt and S used here, no 
crossover to crazing is observed. 
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I. INTRODUCTION 



Adhesion at polymer interfaces is important in many diverse applications such as, col- 
loidal stabilization, [0,0 filler modification of polymeric materials and lubrication 0] and 
enhancing mechanical properties of polymer blends. The interface of interest can be 
either: (i) between two homopolymer melts (A+B) or (ii) between a homopolymer melt 
and a hard surface (A+substrate). Most polymer blends do not mix, because even a weak 
repulsive interaction causes the system to phase separate as the chain length increases. This 
phase separation reduces entanglements at the interface and the adhesive strength is then 
due to relatively weak van der Waals forces alone. For case (ii), the interfacial strength 
depends on the number of chemical bonds between monomers of the polymer and sites on 
the substrate. In the absence of strong chemical bonding, the adhesive strength is again 
dominated by weak van der Waals attractions. For both cases certain additives or compati- 
bilizers can increase the adhesive strength of the interfaces. For case (i), the additive can be 
an AB diblock copolymer. The A(B) block of the copolymer can form entanglements with 
the A(B) melt side respectively. In case (ii), the additive is a chain of the A species with a 
functionalized end-group able to react with the substrate, thereby forming a strong chemical 
bond. The tethered chain of the A species can then become entangled with the melt. In 
both cases, larger numbers of polymer entanglements result in increased interfacial strength. 
The systems simulated in this paper consist of substrate-tethered chains in contact with a 
polymer melt (ie case (ii)), but the results can easily be generalized to polymer/polymer 
interfaces. 

Two key parameters govern the amount of interfacial entanglement, the length (A^t) 
and the areal density (S) of the tethered chains. Adhesion enhancement due to a tethered 
polymer layer shows a surprising non-monotonic behavior as a function of At and S due to 
the phase behavior for the tethered chain layer. The increased work of adhesion W* 

for an elastomeric system versus areal density of graphed chains S is shown schematically 
in Fig. ^ W* is defined as the extra work needed to separate two surfaces in addition to 
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the work exerted against dispersion forces alone. In Region I where each tethered chain 
is isolated from its neighbors, the work of adhesion increases as E increases. However as 
S continues to increase, the tethered chains begin to overlap and phase separate from the 
melt (region II). In this regime, the interpenetration of the tethered chains into the polymer 
melt decreases with increasing E, and as a result W* decreases. At sufficiently high areal 
densities, the polymer melt is completely expelled from the tethered chains thereby causing 
the work of adhesion to fall off to the bare value due to dispersion forces only (region III). 
W* has a similar dependence on Nt. This non-mo notonic behavior has been observed by 
Leger and co-workers P,p^ in elastomeric materials. Kramer and co-workers ||T0|-[14| have 
studied the effects of tethered chains in glassy polymer systems where the effects of crazing 
could complicate the dependence of W* on S. While Fig. |l| relates the macroscopic work 
of adhesion to tethered-chain parameters, it does not address the microscopic processes 
involved in adhesion failure at the melt/tethered-chain interface. This paper investigates 
the effects of Nt, S, T, kg and v on these failure mechanisms. 

The interplay between the microscopic failure mechanisms of tethered chain pullout and 
scission (including their possible relationship to crazing) are not fully understood, partly due 
to the difficulty of direct experimental observation of these phenomena. Molecular dynamics 
(MD) simulations of fracture in highly crosslinked systems ||T6| , p!7|| and crazing ||T8| , p!9| have 
helped to elucidate the crossover from adhesive to cohesive failure of polymer adhesives near 
walls without end-tethered chains. Simulations of tethered chains on small, highly simplified 
models in 2D have been performed, |2O|,0 but were unable to study the effects of chain 



scission in particular. This paper extends the results on our earlier work p2[ of large-scale 
simulations to study the adhesive failure mechanisms of end-tethered chains in contact with 
an entangled polymer melt. Section II contains details related to the MD simulations and 
the method used to implement breakable bonds. In Sec. Ill results of the MD simulations 
are presented and conclusions are given in Sec. IV. 
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II. MD MODEL/METHOD 



We perform continuous-space, molecular dynamics (MD) simulations on a coarse-grained 
model of polymer chains. The polymers are represented by attaching A^t spherical beads of 
mass m with breakable springs. Bead trajectories are obtained by stepwise integration of 
Newton's equations of motion (EOM) using a velo city- Ver let algorithm with a time 
step At=0.006r, where r=cr(m/e)^/^, with a and e setting the length and energy scales 
respectively. The EOM includes terms for a weak stochastic force and a viscous drag force 
with a coefficient on the viscous force term of 0.5r^^. [25| The addition of these two forces 
to the EOM effectively couples the system to a heat bath. Each pair of monomers separated 
by a distance r interacts through a standard (12-6) Lennard- Jones potential Ui^j{r), 



CLj(r) 
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with rc=2.2cr. For simplicity, the substrate interaction is modeled as a flat wall by an 



integrated LJ potential 
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with 2^^"=2.2cr, ewaii=2.0e (strongly attractive wall, no tethered chains) and ewaii=0.1e 
(weakly attractive wall, tethered chains). 

To study the effect of chain scission on adhesion, the standard finite extensible non-linear 
elastic (FENE) potential is altered to allow for broken bonds along a polymer chain. The 
total non-breakable potential Unh{r) (including the FENE potential), between two adjacent 
beads on the same chain separated by a distance r takes the form 

I U,jir) + -0.5KKi\n 1- (^ '^ 

f/nb(r) 



oo 



r < Ro 
r>Ro 



while the breakable potential Uh{r), takes the form 
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UJr) 



Uhjir) + kr^ [(r - ri)(r - r2)] + Uc r < rbr 
r > rbr 

with rbr the length at which a chain bond is defined as broken. MD simulations on a 



similar coarse-grained model p4[ using the FENE potential provide details of the chain 
conformations needed to construct the appropriate starting states. The parameters in the 
bond-breaking potential come from fitting U]j{r) to the region of the aforementioned FENE 
potential {K=30.0e/a'^ and Ro=l.25(x) near its mean equilibrium bond length. This fitting 
procedure results in rbr=l-21cr and a bond-breaking barrier Af/b ~ 20e. The remaining pa- 
rameters in Ub{r) are A:=-409.12e/cr^ ri=1.2a, r2=1.219(T and f/c=42.059e. Fitting Ub{r) 
to a FENE potential this way allows constructing an initial configuration with chains of ap- 
proximately the correct radius of gyration Rg by utilizing the aforementioned chain-building 
algorithm. p4| Figure ^ compares the standard non-breakable potential to the breakable 
potential used in this study. The form of Uh{r) allows for two extrema in the bond poten- 
tial; one stable, global minimum near the FENE potential energy minimum, and one local 
maximum near rbr where the bond force becomes zero. This allows bonds to be removed 
safely from the force calculation without causing large recoil velocities on the resulting chain 
ends. 

The model outlined above results in completely fiexible chains. Experimentally, real 
polymer chains are not completely fiexible but have some degree of stiffness. Chain stiffness 
is added to the model by adding a bending potential energy for each triplet of consecutive 



particles along a chain. The bending potential takes the form |25 

Ue{e) = ke{l + cose) 

with an equilibrium angle of 6o=7t where 6 is defined in Fig. ^. When a bond between two 
monomers is broken during a simulation, all of the triplets associated with that bond are 
removed from the force calculation. 

The initial chain configurations must be constructed appropriately for the model and 
potential parameters chosen. For the systems of very long, entangled chains used in this 



study, it is not feasible to equilibrate the chains by brute-force MD simulation. In gen- 
eral, it is not known how to construct the correct configurations for end-tethered chains 
next to a long-chain melt. To avoid difficulties in initializing systems for the various wet 
brush regimes we restrict the present study to areal densities S in the so-called "mushroom 
regime" (i.e. Region I in Fig. |T]). In this regime each tethered chain interacts weakly 
with other tethered chains and may be constructed as a Gaussian chain. So, for a given Nt 
and S, all chains are constructed as random walks with the correct mean value of Rg. The 
tethered chains are attached to the substrate wall in random locations. The correct value 
of Rg depends on kg, which is known from recent bulk simulations. ||26| The system size is 
adjusted so that the tethered chains do not bridge the box, thereby interacting with the 
opposite wall. A soft potential is used initially to remove overlaps prior to switching on the 
full LJ potential between monomers. The size of the simulation cell is adjusted until the 
pressure P ^ resulting in an overall monomer density of p ~ 0.85o"~^ (p ~ 0.88o"~^) for 
the highest (lowest) temperatures used. For all simulations, the number of beads in each 
of the melt chains is A^m=2500 and the number of tethered chains is nt=30. The chain 
stiffness is taken to be either kg=0 (fully flexible), or k0=1.5e for which the entanglement 
length as calculated from the plateau modulus is Ne=72 and 27 monomer units respectively. 



27| , |28[| The temperatures used range from T=1.0e//cB, which is well above the glass transi- 



tion temperature for fully flexible chains Tg=0.5 — O.Qe/k^, |]T8[ down to T=0.3e/k^. Areal 
densities as low as S=0.002cr^^ with T=1.0e/kB and f =0.0167o"r^^ were used without any 
significant differences with runs using the same parameters for S=0.008(7^^. These results 
confirm the estimates of S which place a system in the mushroom regime. Therefore, to 
increase computational speed the majority of the simulations (all the data presented in this 
paper) use E=0.008o"~^. For systems with crosslinking, the potential for the crosslinks is 
the same as the interactions between monomers bonded along a chain. The crosslinks are 
randomly distributed in the system at an average distance along each chain of approxi- 
mately 70 monomers, which is comparable to Nq. The crosslinks are added only to the melt 
chains after the chain configurations have been constructed and equilibrated. Simulation 



are performed using the massively parallel MD code LAMMPS (suitably adapted to 
include chain scission) developed at Sandia and run on the ASCI Red Terafiop machine and 
Computational Plant (Cplant) clusters. 

Figure ^ shows chain configurations at different times from a tensile pull simulation con- 
sisting of approximately 7 x 10^ particles with A^t=100. Figure ^ shows chain configurations 
at different times from a tensile pull simulation consisting of approximately 2x10^ particles 
with A^t=250. Our largest systems contain close to 10^ particles. The tensile pull is achieved 
by moving only the bottom wall at constant velocity. The tethered chains are bonded to the 
bottom wall of the simulation cell and the top wall has no tethered chains. The interaction 
strength of the top wall is set to be sufficiently strong so that no adhesive failure occurs on 
the top wall during pulling. The interaction of the bottom wall with the melt has a very 
weak attractive component, so that the adhesion enhancement due to the tethered chains 
may be studied independently of the adhesion to the bare wall. The z axis is normal to both 
walls and periodic boundary conditions are used in the x and y directions. The size of the 
simulation cell in the z direction is set to be approximately three times the radius of gyration 
of the tethered chains prior to pulling. The red chains are tethered, blue chains represent the 
melt and green chains were initially tethered and have broken sometime during the tensile 
pull. Figures § and | show the qualitative dependence of the tethered chain dynamics on 
A^f For small At? almost all of the tethered chains are completely pulled out of the melt. 
For relatively large At, all of the tethered chains break somewhere along their length before 
being totally extracted from the melt. 



III. RESULTS 

During the simulation we measure the total work exerted and the number of bonds 
broken by pulling the bottom wall. If wif) is the instantaneous work at time t then the total 
integrated work at t is 

Jo 



To quantify the degree of scission the remaining length of the tethered chains is averaged at 
the end of a pull simulation. The average fractional length of the remaining chains is 

1 n/- 

(F) = - V 131 

where is defined to be length of the ith tethered chain at the end of a pull simulation 
and rit is the number of chains tethered to the pulling wall. The value of (F) quantifies the 
amount of scission in the tethered chains at the end of a pull. {F)=l corresponds to pure 
chain puUout and {F)=0 to chain scission at the wall for each chain. 

A. Chain length dependence 

Figure ^ shows the integrated work over time as a function of time W(t), at different 
values of A^t for T=0.3e/A;B, which is well below Tg. For each tethered chain length, W(t) 
is measured for a sufficiently long time such that each of the tethered chains have either 
completely pulled out or broken. The plateau in Wit) at large t signifies the complete 
debonding of the pulling surface from the entangled melt. Surfaces with longer tethered 
chains require more work to completely pull away from the polymer melt than do short 
chains. However, the maximum work required for surface debonding saturates for large N^. 
This effect can be explained by examining the chain length dependence on tethered chain 
scission. 

Figure |^ shows (F) versus A^t for two values of f at a temperature well above the 
glass transition. For A^t ~ 50 the chains completely pull out of the melt with no breaking. 
For Nt ^ 50 only a fraction of each tethered chain pulls out of the melt before breaking. 
For the longer tethered chain lengths the value of (F) first decreases and then saturates, 
consistent with the data for W(t) shown in Fig. ^ This data shows the entanglement length 
of the tethered chains to be an important length scale. For the fully fiexible model (A;6»=0.0), 
Ne = 72 which is consistent with a value between 50 and 100 for the location of the crossover 
to chain scission as seen in Fig. |^. Figure ^ shows (F) versus A^t for two values of v for 
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T well below Tg. This data has the same qualitative behavior on Nt as the data in Fig. |^a, 
however it is less dependent on the pulling velocity v and the crossover to scission occurs at 
a tethered chain length close to A^e for both values of v. The normal mode relaxation time of 



a chain in an entangled polymer melt increases as the wavelength of the mode increases, pi 
Therefore as the pulling speed v becomes slower, a larger average fraction of each tethered 
chain is able to pull out of the melt before breaking. As T is lowered below the glass 
transition however, the dynamics are dominated by the increased monomeric friction and 
the velocity dependence on chain scission is considerably weaker. 

This is more clearly illustrated in Fig. ^ which shows (F) dependence over two decades 
in velocity T above and below Tg. For T < Tg, the data is nearly independent of the pulling 
speed. For T > Tg, Fig. ^ includes data for an uncrosslinked melt as well as a crosshnked 
network. In most applications, polymer adhesives are crosslinked when T is above Tg in 
order to improve mechanical properties. For large v, the data show that the amount of 
chain scission depends little on whether or not the polymer melt is crosslinked because the 
melt chains cannot relax over short times. This weak dependence of (F) on crosslinking at 
large v can also be seen in Fig. ^ (open squares). Presumably, if v is decreased further 
in the uncrosslinked system for T > Tg then the chains would puUout completely of the 
polymer liquid. However for the crosslinked network at T > Tg, the value of (F) should 
begin to saturate for the lowest pulling velocity, albeit with less chain scission than for all 



the systems below Tg. Extensive simulations on entangled melts have determined that 
the entanglement time t^, ^ 1400r, which is the average time needed for a chain to crossover 
from Rouse-like to reptation-like dynamics. The smallest values of v in Fig. ^ correspond to 
pullout times (i.e. the time required to completely separate the pulling wall from the melt) 
that are order 100 times larger than Tg. 

To further investigate the role of entanglements on chain dynamics, is reduced by 
adding a finite bending energy to the chains. Figures ^ and show (F) versus A't for 
two values of kg respectively, at temperatures above/below Tg. This data has the same 
qualitative behavior on A^t as the data in Figs. |^ and however the value of A^t where 
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the onset of chain scission occurs is smaller for chains with non-zero bending interactions. 
This is consistent with the fact the A^^e is reduced as chain stiffness is increased. The data 
in Fig. 1^ and are from simulations using v=0.167aT^^ which is the faster of the two 
pulling velocities seen in Figs. 0. However, Fig. |^ shows that for T < Tg the chain scission 
data is nearly independent of v. Even for T > Tg, the effect of finite chain stiffness on Ng 
(and its effect on the crossover to chain scission) is clear at the higher value of v. 



B. Temperature dependence 



Figure |T0| shows the integrated work versus time VV(t), as a function of T for A^t=100 
with ke=0.0. For this tethered chain length, the system is only beginning to enter the chain 
scission regime for all values of T. As for Fig. ^, the data are shown out to sufficiently long 
times to show a plateau in VV(t). The data show a steady increase in the plateau value of 
VV(t) as T increases, contrary to what one might expect as T is lowered below Tg. Since 
almost no breaking occurs for this value of N^, the increase in work with decreasing T can 
be attributed solely to an increase in the monomeric friction coefficient. 

Figures [Tl]a and |Tll b show the integrated work versus time W(t), as a function of T for 
A''t=250 with kg=0.0 and 1.5e respectively. For A^t=250, the system is well into the chain 
scission regime for all values of T. The work performed, at early times (i.e. t ^ 2500r) 
increases monotonically as T decreases for both values of kg. However, the plateau values 
of yV{t) show non-monotonic behavior with decreasing T due to scission in the tethered 
chains. Chain scission occurs throughout the pull simulation, from early times oc t^, to 
times corresponding to the beginning of the VV(t) plateau. The smaller the time necessary 
to allow for every tethered chain to have broken, the sooner the plateau in W(t) occurs. The 
non-monotonic behavior, as well as the large jumps, in the plateau values for W(t) suggest 
these effects are not due to changes exclusively in the monomeric friction coefficient since 



this behavior is not present for tethered chains that do not break (see Fig. [10). The plateau 



values for T=0.6-1.0e//cB are nearly identical with a large gap between them and the higher 
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plateau values for all T < 0.6e/k-B- For /ce=1.5e, the plateau value in W(t) has a large jump 
between T=1.0 and O.Ge/k^. However, for T < O.Ge/k^, the plateau values cluster at a lower 
value of VV(t). For A^t > ^e, the effects of the increasing monomeric friction with decreasing 
T in addition to the longer relaxation times for entangled chains combine to produce the 
unexpected, non-monotonic temperature dependence in the W(t) plateau values. 

The temperature dependence of chain scission is shown explicitly in Fig. [l^, which shows 
(F) versus T for the systems shown in Fig. |10| and Figs. |ll]. When A^t=100, negligible chain 
scission occurs for all T. For A^t=250 and T > Tg, the amount of chain scission is constant 
(within errors) for both values of kg. As T is lowered for T < Tg, (F) decreases, which 



helps explain the behavior W(t) in Fig. |Tl|. More work is done because of the increased 
monomeric friction at low T, however the increase in chain scission places an upper limit on 
the total amount of work that can be done pulling the tethered chains before they break. 
For all T, chains with finite stiffness exhibit more scission than fully flexible chains which is 
consistent with the data for (F) in Fig. | which support the fact that is reduced as the 
bending energy increases. 

In addition to monitoring the amount of chain scission in the tethered chains, the number 
of broken bonds on melt chains has also been measured. For the fully flexible chains, kg=0.0, 
a large number of bonds in the melt are broken during a pull simulation for T > 0.6e/k-B- 
However, for T < 0.6e//cB the number of melt chain bonds which break drops suddenly 
to zero. For k0=1.5e, a large number of bonds in the melt are broken for T=1.0e/kB, the 
number dropping to zero for T < l.Oe/k^. The explanation for this behavior is unclear, 
however there seems to be a correlation in the work done and breaking in the melt chains. 
Specifically, the large jumps in the plateau values of W(t) (Fig. |ll]a and pJ]b) occur when 
the number of broken melt bonds drops suddenly to zero. 
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IV. CONCLUSION 



In this study we present results from the first large-scale MD simulations to study the 
dynamics of tethered chains in a dense polymer melt and their effects on adhesion in a 3D, 
realistic polymer model. Data is presented which shows a crossover from pure chain pullout 
to chain scission which depends on the length of the tethered chains relative to the bulk 
entanglement length. For T < Tg, the length of the tethered chains near this crossover to 
scission is consistent with (with/without bending stiffness) for the bead-spring model 
used in this study. This result agrees with experiments on reinforcing the interface between 
two immiscible, glassy polymers with block copolymers. These experiments show evidence 
of copolymer scission and large increases in the work of adhesion, when the molecular weight 
of the copolymer is comparable to the entanglement length. 

However, the experimental data also suggest that in some cases the increase in the work 
of adhesion is due to crazing. The implication is that the copolymer block entangled in 
the melt transfers sufficient stress to initiate crazing. But crazing is not observed in our 
simulations with small values of the wall interaction energy Cwaii even for large Nt and 
T < Tg. Simulations using small values of ewaii correspond to experiments on immiscible 
glassy polymers, because in both cases the attractive forces across interface are too small 
to initiate crazing in the absence of compatibilizers. Simulations were also performed with 
larger values of Cwaii to investigate the inverse problem: could the presence of crazing (due to 
a strong interaction between the melt and the wall) alter chain scission? These simulations 
showed that chain scission is virtually unaffected by the presence of crazing (ie voids) in the 
melt. The lack of crazing with small e^aii might be because of the relatively low tethered 
chain areal densities simulated. Crazing could occur at small ewaii with values of S larger 
than those in the mushroom regime but still low enough to stay out of the completely phase 
separated regime (see Fig. |I]). 

Future work will include developing algorithms to efficiently equilibrate systems with 
tethered chains that are partially phase separated to study the possible effects on crazing. 
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Also, these simulations are being extended to include realistic potentials so that adhesion 
and crazing may be studied in real materials. 

Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Mar- 
tin Company, for the United States Department of Energy under Contract DE-AC04- 
94AL85000. 
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FIGURES 

FIG. 1. Schematic of the increased work of adhesion W* for an elastomeric system versus 
areal density of tethered chains S. (Region I) Mushroom regime, (Region II) Partially overlapped 
(wet brush) regime, (Region III) Phase separated regime. W* shows a similar dependence on Nt. 

FIG. 2. Comparison of the breaking potential Uy^ir) and the unbreakable FENE potential 

FIG. 3. Schematic for bending potential. The three labeled monomers define the triplet used 
in calculating the bending force associated with the angle 9. 

FIG. 4. (color) Chain configurations at two times during a tensile pull simulation for A't=100 
with T=0.3e/fcBi puH velocity v=0.1Q7aT^^ and kQ=0.0. Elapsed times shown are (a) SOOr and 
(b) 600r. The red monomers belong to tethered chains and blue monomers belong to melt chains. 
Green monomers belong to sections of tethered chains which are no longer attached to the bottom 



substrate. Raster3D ]3^ is used to render the images. 



FIG. 5. (color) Chain configurations at three times during a tensile pull simulation for A't=250 
with T=0.3e/fcB) pull velocity v=0.1Q7aT^^ and k0=O.O. Elapsed times shown are (a) 300t, (b) 
900t and (c) 1200t. Coloring scheme is the same as for Fig. 

FIG. 6. Total integrated work W(t) vs. time at different tethered chain lengths A't for 
r=0.3e/A;B, f =0.0167c7t-i and ke=0.0. 

FIG. 7. Average fractional length of remaining tethered chains (F) vs. initial tethered chain 
length A't for (a) T=1.0e/kB, (b) T=0.3e/A;B- The symbols denote v=0.l67aT~^ (filled square) and 
t'=0.0167(7T^^ (filled triangle). The open squares show similar results for a crosslinked network 
with ~ 18 crosslinks per melt chain for u=0.0167o"t^^. For all data kg=0.0. 
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FIG. 8. Average fractional length of remaining tethered chains (F) vs. tensile pull velocity v 
for Ni=250 at T=1.0e/kB for the crosslinked (open triangles) and uncrosslinked (solid triangles) 
system and at T=0.3e//cB for the uncrosslinked system (solid squares). For all data fce=0.0. 

FIG. 9. Average fractional length of remaining tethered chains {F) vs. initial tethered chain 
length Nt for (a) T=1.0e/kB, (b) r=0.3e/A;B- The symbols denote k0=O.Q (filled square) and 
ke=1.5e (filled triangle). For all data v=0.167aT~^. 

FIG. 10. Total integrated work W{t) vs. time as a function of temperature for A?t=100, 
w=0.0167c7T-i and ke=0.0. 

FIG. 11. Total integrated work W(t) vs. time as a function of temperature for Ni=250 and 
u=0.0167crr-i. (a) k0=O.O (b) ke=l.be. 

FIG. 12. Average fractional length of remaining tethered chains {F) vs. temperature T for 
i;=0.0167crr-^ The symbols denote Nt=lQQ, ke=Q.Q (filled square), Art=250, ke=Q.Q (filled trian- 
gle) and A^t=250, k0=l.he (open triangle). 
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